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Abstract. Spin-dependent electronic transport through a quantum dot has been analyzed theoretically 
in the cotunneling regime by means of the second-order perturbation theory. The system is described by 
the impurity Anderson Hamiltonian with arbitrary Coulomb correlation parameter U. It is assumed that 
the dot level is intrinsically spin-split due to an effective molecular field exerted by a magnetic substrate. 
The dot is coupled to two ferromagnetic leads whose magnetic moments are noncollinear. The angular 
dependence of electric current, tunnel magnetoresistance, and differential conductance are presented and 
discussed. The evolution of a cotunneling gap with the angle between magnetic moments and with the 
splitting of the dot level is also demonstrated. 

PACS. 72.25.Mk Spin transport through interfaces - 73.63.Kv Quantum dots - 85.75.-d Magnetoelec- 
tronics; spintronics: devices exploiting spin polarized transport or integrated magnetic fields - 73.23.Hk 
Coulomb blockade; single-electron tunneling 



1 Introduction 



The vast step forward in fabrication of tunnel junctions 
has enabled the implementation of extremely small metal- 
lic grains and semiconductor quantum dots coupled through 
tunnel barriers to external reservoirs |1I2) . Currently, trans- 
port through such ultra-small devices is being extensively 
studied because of many future application possibilities 
and, more importantly, because of beautiful physics that 
emerges in these nanoscale systems. In nanostructures not 
only the manipulation of a single electron charge is pos- 
sible, but - when coupled to ferromagnetic leads - also 
the manipulation of a single electron spin. This is why 
those systems are considered to play an important role in 
spintronic devices. 

An interesting feature of electronic transport through 
nanoscale systems coupled to ferromagnetic leads is the 
tunnel magnetoresistance (TMR). The TMR effect, al- 
though discovered already three decades ago in planar 
junctions [5], is still of current interest. It consists in a 
change of the system conductance when relative orienta- 
tion of the magnetic moments of external leads switches 
from antiparallel to parallel alignment. In a general case, 
magnetic moments of the two electrodes can form an ar- 
bitrary angle (noncollinear configuration). Qualitatively 
different features of electron transport through nanosys- 
tems appear due to discrete charging effects, leading to 
Coulomb blockade and Coulomb oscillations of electric 
current. The interplay of charge and spin effects gives rise 



to new interesting phenomena, like for instance TMR os- 
cillations with applied bias and gate voltages [41516171813] . 

Theoretical considerations of electron tunneling through 
quantum dots attached to ferromagnetic leads have al- 
ready been reported in a number of papers. Most of them, 
however, are restricted to spin-dependent transport in sys- 
tems with collinear alignment of the electrodes' magneti- 
zations. Basically, all transport regimes in such a geometry 
have already been addressed, including sequential (first 
order) tunneling |10lll| . cotunneling (second order) [T2]. 
resonant tunneling |13I14| . and Kondo regimes |15I16| . 
Spin polarized transport through quantum dots coupled to 
ferromagnetic electrodes with noncollinear magnetic mo- 
ments is still not fully explored, although it has already 
been considered in a couple of papers |17I18I19I2UI21I22| . 

In this paper we address the problem of second-order 
(cotunneling) spin-dependent transport through quantum 
dots coupled to ferromagnetic leads with arbitrary con- 
figuration of the in-plane magnetic moments of external 
electrodes. The considerations are limited to single-level 
quantum dots. Moreover, the level of the corresponding 
isolated dot is assumed to be spin-split due to coupling 
between the dot and a magnetic substrate on which the 
dot is deposited. The corresponding splitting is assumed 
to be larger than the level width r due to coupling of the 
dot to external leads. 

The paper is organized as follows. The model and method 
are described in section II. Transport through an empty 
quantum dot in the cotunneling regime is described in 
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section III, where also the relevant numerical results are 
presented and discussed. Section IV covers the problem of 
cotunneling through a singly occupied dot. Final conclu- 
sions are given in section V. 



2 Description of the model and method 

The system considered in this paper consists of a single- 
level quantum dot coupled through tunnel barriers to two 
external ferromagnetic leads, whose magnetizations are 
oriented arbitrarily in the plane of the structure. The dot 
is assumed to be deposited on a ferromagnetic substrate 
which strongly interacts with the dot and leads to spin- 
splitting of the dot level. The splitting is assumed to be 
larger than the level splitting due to exchange interaction 
between the dot and electrodes. The exchange interaction 
results from tunneling processes and is of the first order 
in the coupling parameter r |18I19| . When neglecting the 
exchange coupling between the dot and leads, one may as- 
sume that the level splitting is constant, i.e., independent 
of applied voltage. 

Coupling of the dot to ferromagnetic substrate is de- 
scribed by an effective molecular field B s , lying in the 
plane of the structure. Thus, the magnetic moments of 
the electrodes and the molecular field are all in a com- 
mon plane. The effective molecular field determines the 
global quantization axis for electron spin on the dot (axis 
z antiparallel to the molecular field). In turn, the local 
quantization axis in the leads is determined by the corre- 
sponding local magnetization direction, and is parallel to 
the net spin of the lead (thus, being antiparallel to the lo- 
cal magnetic moment). Furthermore, we assume that the 
net spin of the left (right) lead forms an angle ipi, (<^r) 
with the global quantization axis as shown in Fig. ^ In 
order to distinguish between different quantization direc- 
tions, the majority (minority) spins in the local reference 
systems are labelled with a = +(— ), while spin projection 
on the global quantization axis is denoted as a =] (J.) for 
spin- up (spin-down) electrons. 

Hamiltonian of the system has the general form H = 
Hl + Hn + Hd + Ht, where Hl and describe the 
left and right leads as reservoirs of noninteracting quasi- 
particles, Hy> is the dot Hamiltonian, and tunneling pro- 
cesses between the electrodes and dot are included in Ht- 
The lead Hamiltonians are diagonal in the respective local 
coordinate systems, H a = J2cr=+- J2kea <^akaO i akci a aka 
(for a — L, R), with e a ka being the energy of a single elec- 
tron with wavevector k and spin a in the lead a, whereas 
a aka an d a aka denote the corresponding creation and an- 
nihilation operators. 

The dot is described by the Anderson Hamiltonian, 
which in the global reference frame can be expressed as 
i?D = So-=T l £ <y^a^<y + Un^n^, where e a is the energy of 
an electron with spin a, e a = e=F <?/iB-B s , and d\, (d a ) cre- 
ates (annihilates) a spin-cr electron. Here, g is the Lande 
factor and e is the dot level energy in the absence of 
magnetic electrodes and molecular field due to the sub- 
strate. The second term of the dot Hamiltonian describes 
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Fig. 1. The schematic of a magnetic quantum dot coupled to 
ferromagnetic leads. Classical spin of the lead a (a = L, R) can 
form an arbitrary angle ip a with the dot spin quantization axis 
(axis z), as indicated. 



Coulomb interaction of two electrons of opposite spins re- 
siding on the dot, with U denoting the corresponding cor- 
relation energy. The tunnel Hamiltonian takes the form 



+ {T ak+ al k+ sin -j- + T ak ^a f ak _ cos -y)d; + h.c.p) 

with T akt7 denoting the tunnel matrix elements between 
the dot states and majority (o = +) or minority (p = — ) 
electron states in the lead a when (p a = 0. 

Due to the coupling between dot and leads, the dot 
level acquires a finite width. When the magnetic moment 
of lead a and the molecular field acting on the dot are 
parallel, the corresponding contribution r° to the dot 
level width may be written as r° = 2ir\T aa -\ 2 p acr , where 
p aa is the spin-dependent density of states for the ma- 
jority (er = +) and minority (a — — ) electrons in the 
lead a. The parameters will be used in the following 
to parameterize strength of the coupling between the dot 
and lead a. It is convenient to express the coupling pa- 
rameters in terms of spin polarization defined as P a = 
( r a - r~)/(r+ + r~). Thus, the coupling strength can 
be written as F± = r a (l±P a ), where r a = {r++T-)/2. 
In our considerations we assume J\, = -Tr = r/2. 

In the case of ferromagnetic leads, the coupling of 
the spin-up dot level is different from the coupling of 
the spin-down level, which is due to different densities of 
states for spin-majority and spin-minority electron bands 
in the leads. This may result in the splitting of the dot 
level (12123) . Here, we assume that the dot level splitting, 
A = si — £f, due to the molecular field is larger than the 
coupling parameters, A 3> -T^. In other words, we assume 
that the exchange interaction between the dot and mag- 
netic leads is much smaller than the Zeeman energy due 
to the molecular field and can be neglected. An electron 
residing on the dot has then either spin up or down. Thus, 
the corresponding density matrix in the global quantiza- 
tion system is diagonal in the spin space In that case, 
only the second-order processes (with respect to the tun- 
neling Hamiltonian) have to be taken into account in the 
Coulomb blockade regime. When the above condition is 
not fulfilled, the exchange interaction cannot be ignored 
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and therefore the first-order processes, which are respon- 
sible for the exchange coupling (and do not contribute to 
charge transport), have to be considered. 

In the Coulomb blockade regime, the energy needed to 
put an electron on the dot surpasses the energy provided 
by the transport voltage and the sequential tunneling is 
exponentially suppressed. The dot is then in a well de- 
fined charge state, and quantum charge fluctuations are 
suppressed. Although the energy conservation prohibits 
the first-order tunneling transitions, the current can still 
be mediated by higher-order tunneling processes involving 
correlated tunneling of two (cotunncling) or more elec- 
trons via intermediate virtual states [251 ■ 

The rate of electron cotunneling from a spin-majority 
state in the lead a to a spin-majority state in the lead 
(3 can be determined using the second-order perturbation 
theory |25I26| . and is given by 



'a/3 



2vr 

T 



E 



($+\H T \<l> v ){cl> v \H T \$p 



6(Ei-Ef), 



(2) 

where |^+) and \@t) are the initial and final states of the 
system, whereas \<P V ) is a virtual state. The corresponding 
energies are denoted as £,-, £/, and e v . 

One can distinguish cotunneling processes that change 
the magnetic (and consequently also energetic) state of 
the dot, and processes affecting neither magnetization nor 
energy of the dot. The former (latter) processes will be re- 
ferred to as inelastic (elastic) ones. The elastic cotunneling 
processes are fully coherent and do not change the dot 
occupation probabilities. Contrarily, inelastic cotunneling 
influences the occupation numbers of the dot, and can take 
place only when the dot is occupied by a single electron 
- either spin-up or spin-down. Furthermore, one can also 
distinguish between single-barrier and double-barrier co- 
tunneling. Only the latter processes contribute directly to 
the current. However, the inelastic single-barrier cotunnel- 
ing processes can change the occupation probabilities, and 
consequently can also influence the electric current flowing 
through the system. In particular, inelastic single-barrier 
cotunneling processes which reverse spin of the dot play 
a significant role. This is because they can open system 
for more efficient tunneling processes, when the system is 
blocked by an electron of a given spin orientation residing 
on the dot |T2] . 

In the following we will consider two different situa- 
tions. The first one corresponds to an empty dot (due to 
particle-hole symmetry the results can be adapted to the 
case of doubly occupied dot). The second situation is the 
case of a singly occupied dot. 



3 Cotunneling through an empty dot 

When the dot level is far above the Fermi energy of the 
leads, Sfj ^> k&T, r, \eV\, there are no electrons on the 
quantum dot and electric current can flow only due to 
elastic cotunneling processes. Because of the particle-hole 
symmetry, a similar analysis can be directly performed for 



a doubly occupied dot, when e a + U «C and \E a + U\ 3> 
k^T, r, \eV\. As in the case of an empty dot, the current 
can then flow only due to elastic cotunncling. 



3.1 Theoretical description 

Electric current / flowing from the left to right lead is 
given by 



E fa 



7lr,o — 7rl,o 



<7,£T'— + , 



(3) 



with 7lr*x b em g the elastic cotunneling rate for transi- 
tion from the left to right leads, when the dot is in the 
state \$ x ) [x = in Eq.(3)] and when a majority (cr = +) 
or minority (a = — ) electron of the left lead tunnels to 
majority (a' = +) or minority (cr' = — ) electron band in 
the right lead. Similarly, is the elastic cotunneling 

rate for transition from the right to left electrodes. In Eq. 
(3) — e denotes the electron charge (e > 0). 

The transition rate for electrons tunneling from the 
majority spin band in the left lead to the majority spin 
band in the right lead is given by the formula 



7lr,o 



x|T L . 



2tt 
= Y 
| 2 |T R+ | 2 /(e L )[l 



de L de R pL+Pu- 



/(<*)] 



cos((^ l /2)cos(^r/2) sin(^ L /2)sin(^ R /2) 



£l + PL - £] 
x5(el + PL - £r - Mr) ; 



EL + PL - El 



(4) 



with /il (pr) denoting the electrochemical potential of 
the left (right) lead and /(e) being the Fermi-Dirac dis- 
tribution function, f(e) = l/[exp(e/fcBT) + 1]. We assume 
/il = — eV/2 and /ir = eV/2, and the energy is measured 
from the Fermi level of the leads in equilibrium situation 
(V = 0). The integrals in Eq.(4) can be calculated quite 
easily using the contour integration method, as described 
in the Appendix. Following this procedure one finds 



7 L + rTo + = ^/b(Mr-M L ){ 



cos2 <|w^ 2(£t) 



• 2 fL ■ 2 VR a 1 s 

-sin ysin — A 2 {e l ) 



sin y L sin y? R 

' 2(e T - £i ) 



[AM - Afa)] 



(5) 



where /b(e) is the Bose function, /b(e) = l/[exp(e/fcBT) — 
1], and A n (s a ) = B n (e a - Pr) - B n (e a - Pl), with B n (x) 
defined as 



B n (x) = Re — - $ 



2 27rfc B T 



(G) 



Here \P(z) is the digamma function. Similar formulae can 
also be derived for the other transition rates. 
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Fig. 2. (Color online) The cotunneling current (a,c) and the 
TMR effect (b,d) as a function of the bias voltage (left column) 
and the angle ipn. (right column) . The parameters assumed for 
calculations are: k-gT = P, £-[■ = 18P, £j = 22P, Pl = Pr = 
0.5, V3 L = 0. 



Having found all the cotunneling rates, one can calcu- 
late the electric current, 



e p 2 

Ah 



(1 - P L cos<£l)(1 - Pr coa (p R ) A 2 (e i 



-(1 + P L COS 99 L )(1 + PrCOS (^r)A 2 (£ T ) 

2PlPr sin <y9 L sin </> R 



MeO-Aifa)]}. (7) 



In the following we will present numerical results on 
the electric current and the associated magnetoresistance. 
The TMR effect is described quantitatively by the ratio 



TMR: 



I(<Pl,<Pr) 



(8) 



where /((^l 5 Vr) is the current flowing in the noncollinear 
configuration described by the angles </Jl and tpp, whereas 
Ip is the current flowing in the parallel configuration cor- 
responding to y>L = <pp = 0. 



3.2 Numerical results 

The formula (7) for electric current corresponds to the 
situation, where the magnetic moments of the leads and 
the effective molecular field acting on the dot are oriented 
arbitrarily in the plane of the structure. Further numer- 
ical analysis will be restricted, however, to the following 
two situations: (i) the magnetic moment of the left lead is 
parallel to the molecular field acting on the dot (</?l = 0), 
while the magnetic moment of the right lead can have an 
arbitrary orientation, and (ii) the magnetic moments of 



both leads can rotate symmetrically in the opposite di- 
rections, ifp — —(fh- The case (i) corresponds to the sit- 
uation when magnetic moment of one lead and magnetic 
moment of the layer producing the molecular field acting 
on the dot are fixed along the same direction. This can be 
achieved for instance by a common antiferromagnetic un- 
derlayer with strong exchange anisotropy at the antiferro- 
magnet/ferromagnet interface. The situation (ii), in turn, 
corresponds to the case when both leads are equivalent 
and their magnetic moments can be rotated simultane- 
ously by an external magnetic field. Now, we will analyze 
both situations in more detail and we begin with the case 

(i). 



3.2.1 Case (i): y>L = 

When ifi, = 0, the electric current is given by the formula 

/= ^-[(l-PL)(l-PRCOS^ R )4 2 ( £i ) 

+ (l + P L )(l + P R cos(^ R )A 2 ( £T )], (9) 

which follows directly from Eq. (7). In Figure [3 we show 
the current (a,c) flowing through the system and the cor- 
responding TMR (b,d) as a function of the bias voltage 
for several values of angle ipp, and as a function of the 
angle tpp for several values of the bias voltage. The cur- 
rent decreases and TMR increases as the angle ipp varies 
from ipp = to ipp = 7r, which corresponds to the transi- 
tion from parallel to antiparallel magnetic configurations. 
Parts (c) and (d) show explicitly this angular dependence. 
Both the current and TMR vary monotonously with ipp_ 
(for < <pp < 7r) and electric current reaches minimum, 
while TMR maximum, at ip-R = n, i.e., in the antiparallel 
configuration. Such a behavior is typical of normal spin 
valves and results from spin asymmetry in tunneling pro- 
cesses. It is also worth noting that TMR is only weakly 
dependent on the bias voltage [see Fig. 2(b)]. 

Assuming the same spin polarization of the leads, Pl = 
P R = P, one finds the following explicit formula for TMR 
in the zero bias and zero temperature limits: 



TMR = 



(l-cos^)P (1 + P)e 2 -(1 



P)e 2 



(1 + P)(l +Pcosv>ji)ef + (1 - P)(l - P cos tp R )e% ' 

(10) 

The above formula describes the angular variation of TMR 
and shows explicitly that TMR reaches maximum for ip = 
7r. This maximum value is given by the expression 



TMR (n 



2P 



1-P 2 



P 



(11) 



The first term in the bracket of the above equation gives 
the Julliere's value of TMR, whereas the second term de- 
scribes enhancement of the tunnel magnetoresistance due 
to level splitting. Such an enhancement of TMR may be 
of some interest from the application point of view. 
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local minimum) is given by 
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Fig. 3. (Color online) The cotunneling current (a,c) and the 
TMR effect (b,d) as a function of the bias voltage (left column) 
and as a function of <p = ip-R. = —fL (right column). The other 
parameters are the same as in Fig. 2. 



3.2.2 Case (ii): ipn 



-<PL 



Consider now the situation (ii), when ipn = —ip^ = >P- The 
magnetic moment of the left lead rotates now together 
with the magnetic moment of the right lead, but in the 
opposite direction. The electric current is then given by 
Eq. (7) with ipn = ip and ip-^ = — ¥?■ The corresponding 
bias and angular dependence of the cotunneling current 
and TMR is displayed in Fig. 

It is interesting to note that the angular dependence 
of electric current and TMR differs now from that found 
above for the situation (i). In order to understand this dif- 
ference one should take into account the fact that <^r = 
ipii — corresponds to the parallel configuration, whereas 
the situation with ipn = —tp-L = 7t/2 corresponds to the 
antiparallel configuration with magnetic moments oriented 
perpendicularly to the molecular field acting on the dot. It 
is also worth noting that electric current reaches minimum 
and TMR maximum not exactly at tpn = —tp^ = 7r/2, 
but for the configuration which is close to the antiparallel 
one. Moreover, position of these minima and maxima de- 
pends on the bias voltage, as shown in Fig. 3(c,d). In turn, 
the case </?r, = — (pi, = 7T corresponds again to the paral- 
lel configuration, but with the magnetic moments of the 
electrodes being antiparallel to the molecular field. This 
leads to a local maximum of electric current (minimum in 
TMR) at ipn = —(p-L = 7r. It is interesting to note that the 
two parallel configurations (aligned and anti-aligned with 
respect to the molecular field) are not equivalent, and con- 
sequently the corresponding currents and also TMR values 
are not equal. 

When assuming equal spin polarizations of the leads 
and low bias and temperature limits, one finds that the 
TMR effect at tp-R = — <P~l = tt (which corresponds to its 



TMR ( 



min) 



4P( £ 2 



(1 



(1 



P) 2 £2 ' 



(12) 



In the same limit and for a nonzero spin polarization P, 
one can show that the two maxima of TMR appear at 
ip — arccosa and ip = 2n — arccosa, with a = — A/P(s-\ + 
ejj. In the limit of P = TMR vanishes by definition. 
In the case when \A\ < P(s-\ + £j,), the maxima appear 
approximately at ipn = —iph = 7r/2 and </jr = — = 
3"7r/2. The corresponding TMR value is then equal 



TMR (r 



<) _ P[(2 + P) £x -(2-P) £T ]( £T+£i ) 



2P% T£i 



(13) 



4 Cotunneling through a singly occupied dot 

By applying an external gate voltage to the dot, one can 
tune position of the level energy and this way also the 
dot occupation. When e a is negative and e a + U positive, 
the dot is singly occupied at equilibrium for P, k^T <C 
|£<r|) \sa + U\, and the system is in the Coulomb blockade 
regime. As before, we will consider the two situations (i) 
and (ii) defined in the previous section. 



4.1 Theoretical description 

In the case studied in the preceding section the dot was 
empty, and the second-order current was mediated only by 
elastic cotunneling processes. When the dot is singly occu- 
pied, the cotunneling current can also flow due to inelastic 
cotunneling, in which the electrons tunneling to and off the 
dot have opposite spin orientations (and consequently also 
different energies due to the level spin-splitting). These 
inelastic cotunneling processes determine the occupation 
numbers of the dot. 

The probabilities p-\ and pj , that the dot is occupied ei- 
ther by a spin-up or spin-down electron, can be calculated 
from the following stationary master equation: 



0= E 

a,/3=L,R 



(14) 



and taking into account also the normalization condition, 
Pi +Pi = 1. Here, 7 Q/ 3, (T =^5' denotes the rate of inelastic 
cotunneling from lead a to lead (3 with the simultaneous 
change of the dot spin from a to a = —a. The inelas- 
tic cotunneling processes can take place through one of 
the two virtual states, |#o) an d |^d), which correspond 
to the empty and doubly occupied dot, respectively. The 
corresponding energies are Eq — for the empty state and 
Ed = £1 + E[+U for the doubly occupied state. 

The rate of inelastic processes which transfer a spin- 
majority electron from the left lead to the spin-down level 
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of the dot and a spin-up electron from the dot to the spin- 
majority electron band in the right lead is given by 

TlhvN-I = ~kJJ de L d£ RPL+PR+|7L+| 2 |T R+ | 2 
xsin ^ cos2 ^ /(eL) [l_ /(eR) ] 



1 



1 



x£(£l + Ml - £r - MR + £t ~ £ i)- 



(15) 



The total rate of inelastic cotunneling from the left to 
right leads, which changes the dot state from to 
can be found by summing up over the spin-majority and 
spin-minority electrons, 



7lr,t=h 



E 



7lr,t=H.- 



(16) 



The corresponding analytical expression can be derived in 
a similar way as in the case of elastic cotunneling described 
in the previous section and takes the form 



7lr,t=h 



/b(mr-Ml -£t + £ l 



x 1 jl tT sin — cos h i ■ J n sin — sin — 

V LR 2 2 LR 2 2 



+r T :r R + cos 2 f cos 2 ^ 



2 PL . 2 ¥>R\ 

It r„ cos — sm — i 



L 'R ^ u0 ^ ~~ 2~ L R 

x<l -B 2 (£ T - fin) - B 2 (ej. - Ml) 



-B 2 (£ T + U - /it R ) - B 2 (e; + ET - ml) 
2 

-— [5i(e T -/i R ) -Bi(s| - ml) 

-Bl(£ T + C/ - m ) + -Bi(e; + C/ - ml)] J. 
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(17) 



The cotunneling rate 7lr,|=>t can ^ e calculated in a sim- 
ilar way. In turn, the rate of elastic cotunneling can be 
found as described in section III. 

Having calculated the cotunneling rates and the occu- 
pation probabilities, one can determine the current flowing 
from the left to right leads. The contributions /; nc i and I c \ 
due to inelastic and elastic cotunneling, respectively, are 
given by 



-e ^2 (7LR,cr^a ~ 7RL, C r^a)p C r, 

<T=U 



E 



7LR,<7 



7rl,<t 



<t=Ti1 <j' ,<j" =+■- 
The total cotunneling current I is then equal 

I = hne\ + lei- 



(18) 
(19) 

(20) 



The analytical expression for the current in the case of a 
singly occupied dot is cumbersome and will not be pre- 
sented here. 
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Fig. 4. (Color online) The cotunneling current (a,c) and the 
TMR effect (b,d) as a function of the bias voltage (left column) 
and the angle ifn (the right column). The parameters assumed 
for calculations are: ksT = T, Sf = —22F, ey = — 18F, U = 
40r, P L = P R = 0.5, <^ L = 0. 



4.2 Numerical results 

As far as physics is concerned, the situation with singly 
occupied dot becomes more interesting. As before, we will 
analyze the two situations (i) and (ii). 



4.2.1 Case (i): ip L = 

When the magnetic moment of the left electrode is fixed 
(<Pl = 0) and the magnetic moment of the right lead is free 
to rotate, the angular and bias dependence of electric cur- 
rent and TMR is shown in Fig. 4. Except for the parallel 
configuration, the current-voltage curves are now asym- 
metric with respect to the bias reversal [see Fig. 4(a)]. 
This asymmetry also leads to related asymmetric behav- 
ior of TMR [Fig. 4(b)]. Moreover, for positive bias voltage, 
the TMR effect can change sign and become negative in a 
certain range of the bias and angle values. Such an asym- 
metry in transport characteristics with respect to the bias 
reversal is of some importance for applications, particu- 
larly when the current is significantly suppressed for one 
bias polarization (diode behavior). 

In order to account for the bias asymmetry, let us con- 
sider only the antiparallel configuration. One should then 
realize that owing to the level splitting, the single-barrier 
inelastic cotunneling processes can occur only when the 
dot is occupied by a spin-down electron. This follows sim- 
ply from the energy conservation rule. Thus, the single- 
barrier processes can assist the fastest double-barrier co- 
tunneling processes (spin-up electron tunnels through the 
left barrier and spin-down electron tunnels through the 
right barrier), but only for positive bias. This is because 
the fastest processes can occur when the dot is occupied 
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Fig. 5. (Color online) Differential conductance corresponding 
to the situation shown in Fig. 4(a), calculated for two different 
temperatures. 
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Fig. 6. (Color online) Differential conductance for different 
values of the level splitting for the parallel (a) and antiparallel 
(b) magnetic configurations. The parameters are: k^T = 0.2_T, 
e T = e - A/2, gj. = e + A/2, e = -20r, U = 40r, and 
Pl = Pr = 0.5. 



by a spin-down electron for negative bias and by a spin-up 
electron for positive bias. From this follows that the con- 
ductance is larger for positive than for negative bias volt- 
age. This is indeed the case in the characteristics shown in 
Fig. 4(a). The above described mechanism of the asymme- 
try with respect to the bias reversal does not hold when 
magnetic moments of the two leads are parallel, so the 
corresponding current-voltage curve is symmetrical. 

To demonstrate the above described asymmetry more 
clearly, the corresponding differential conductance is dis- 
played in Fig. 5(a). The asymmetry for tp K > is clearly 
evident. Besides the asymmetry, an additional interest- 
ing feature of the differential conductance is also visible, 
namely the characteristic deep in the small bias regime. 
This deep is a consequence of the suppression of inelas- 
tic double-barrier cotunneling events when \eV\ < \A\. 
For \eV\ > \A\, the inelastic cotunneling processes are al- 
lowed, leading to an enhanced conductance. The asymme- 
try and zero bias anomaly are even more evident at lower 
temperature, as shown in Fig. 5(b). Such a suppression 
of the inelastic cotunneling events at small bias was used 
recently as a spectroscopic tool to determine spin splitting 
of the dot level and the corresponding g-factor |28| . 

The angular variation of electric current and TMR 
reveals further new features. For negative bias there is 
a maximum of absolute value of electric current in the 
parallel configuration and a minimum in the antiparallel 
configuration. For positive bias, however, the electric cur- 
rent has a maximum for noncollinear configuration, as de- 
picted in Fig. 4(c). The nonmonotonous variation of elec- 
tric current with the angle </?r leads to the corresponding 



nonmonotonous variation of TMR, shown in Fig. 4(d). 
Furthermore, TMR may now become negative, as already 
mentioned before. 

Numerical results presented so far were shown for a 
single value of the level splitting, A = 4T. From the ex- 
perimental point of view, variation of the conductance as 
a function of the level splitting (induced for instance by 
a strong external magnetic field) , allows one to determine 
some interesting transport and spectroscopic characteris- 
tics. Therefore, in Fig. 6 we show the differential conduc- 
tance for different splitting of the dot level in the paral- 
lel and antiparallel configurations. By measuring width of 
the conductance deep, one can determine for instance the 
spectroscopic (7-factor [2HI ■ 

In the case of a deep Coulomb blockade regime and 
\eV\,ksT <C |Z\|, one can derive an approximate formula 
for the deep in differential conductance due to the sup- 
pression of inelastic cotunneling. Writing £j = e — A/2 
and £j = e + A/2 one finds then the following expression: 



G = 



r 2 e 2 



(I + P L COS(^ L )(f + P R COS<£r) 



(I-P L COS(^ L )(I 



A/2? 
P R cos v?r) 



{e + U + A/2Y 
8P L P R sin ip h sin ip R 
(e- A/2)(e + U + A/2) 



(21) 



which is valid for arbitrary magnetic configurations. This 
expression approximates the plateaus shown in Fig. 5. 
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When assuming ip^ — <pr — (which corresponds to the 
parallel configuration), Eq. 1)21(1 simplifies to the following 
form 



G P = 



(1 + P L )(1 + P R ) (1-Pl)(1-P r ) 



(e-Zi/2)2 



(£ + <y + z\/2) 2 



(22) 

whereas for the antiparallel alignment (<pi, = 0, <p>R = tt) 
it becomes 



Gap 



r 



2 C 2 



4/7 



(1 + Pl)(1-P r ) , (1-P L )(1 + P R 



{s-A/2f 



(e + U + A/2) 2 



The above two expressions describe the plateaus in dif- 
ferential conductance shown in Fig. 6. It is also worth 
noting that generally Gp > Gap- Furthermore, both Gp 
and Gap vary monotonously with spin polarization of the 
leads - in the case of nonmagnetic leads (Pl = Pr = 0) 
Gp = Gap, whereas for Pl = Pr = 1 (which corresponds 
to half-metallic leads) Gp is maximal and Gap = 0. 



4.2.2 Case (ii): ipR 



Transport characteristics in the second situation, i.e., for 
the case when </?r = — (pi, = are displayed in Fig. 7. One 
can note that the current is now almost independent of the 
magnetic configuration. Nevertheless, the angular depen- 
dence of the current becomes more visible in the corre- 
sponding differential conductance, plotted in Fig. 8(a) for 
different values of the angle (p, and also in Fig. 8(b) for the 
same situation, but for much lower temperature. The co- 
tunneling gap due to suppression of the inelastic processes 
is also clearly visible. The deep in differential conductance 
for \eV\ <C fceP, |^| is given by Eq. l(2"T|) . Since the system 
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Fig. 7. (Color online) The cotunneling current (a,c) and the 
TMR effect (b,d) as a function of the bias voltage (left column) 
and tp = ipn — — <Ph (right column). The other parameters are 
the same as in Fig. 4. 
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Fig. 8. (Color online) Differential conductance for the situa- 
tion shown in Fig. 7(a), calculated for two indicated tempera- 
tures. 



is now symmetric, the current-voltage curves (and conse- 
quently also the differential conductance) are symmetric 
with respect to the bias reversal. 

The TMR effect reaches maximum in the zero bias 
limit, V = 0, as shown in Fig. 7(b). In turn, angular vari- 
ation of TMR reveals two maxima (and also two minima), 
similarly as it was in the case of empty dot, but now the 
maxima appear strictly for tp = tt/2 and ip = 3tt/2. On the 
other hand, one minimum of TMR occurs at tpR = (pi, = 0, 
where TMR vanishes by definition. Tunnel magnetoresis- 
tance vanishes also in the second parallel configuration, 
when both magnetizations are antiparallel to the molecu- 
lar field (cpR — — = tt). This is due to the fact that the 
parameters assumed for numerical calculations correspond 
to a symmetrical Anderson model, i.e., U — — e-f — E[. 
When the system becomes asymmetric, e.g., when U in- 
creases (decreases) while the other parameters are con- 
stant, the minimum in TMR at <p = tt is shifted down 
(up) and when U ^> — £| — TMR has only one maxi- 
mum at ip = it, as shown in Fig. 9(b). On the other hand, 
if U < — £j — £[, the minimum in tunnel magnetoresistance 
at ip — tt becomes negative, which is shown in Fig. 9(b) 
for U — 35P. Thus, by changing the model parameters 
one may significantly enhance or reduce the TMR effect. 

In the limit of \eV\,k B T < \A\ and for P L = Pr = P, 
the TMR ratio for <pn = — = 7T can be expressed as 



TMR: 



4P(2e +U)(U + A) 



(1 - P) 2 (e + U + A/2) 2 + (1 + P) 2 (e - A/2) 2 ' 

(24) 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 



9 




Fig. 9. (Color online) The angular dependence of the differen- 
tial conductance and TMR for several values of the Coulomb 
interaction parameter U at the bias voltage eV = 2F . The 
other parameters are the same as in Fig. 7. 



From the above expression follows that the sign of TMR 
at <p = 7r depends on the ratio e/U. If U = — 2e, TMR 
vanishes, whereas for U ^ — 2e, TMR is positive (nega- 
tive). In the case of U —2e, the TMR effect is given by 
4P/(1-F) 2 . 

The corresponding behavior of the differential conduc- 
tance is shown in Fig. 9(a). Now, the maximum in the con- 
ductance at ip — 7r for a symmetric model changes into a 
minimum for U 3> — £f — ej.. 



5 Conclusions 

We have considered analytically and numerically cotun- 
neling current and associated tunnel magnetoresistance 
through a single-level quantum dot coupled to two ex- 
ternal ferromagnetic leads. The dot level was assumed to 
be spin-split due to an effective molecular field originat- 
ing from a magnetic substrate on which the dot is de- 
posited, and the splitting was assumed to be larger than 
the characteristic parameter _T (level width) describing the 
dot-lead interaction. Although the general formulae were 
derived for arbitrary in-plane orientation of the molecular 
field and of the magnetic moments of the leads, detailed 
numerical analysis was performed for two particular con- 
figurations, which seem to be of particular interest. At this 
point we would like to note, that a related problem has 
been recently studied by Pedersen et al. who con- 
sidered cotunncling current in a similar system, but for 
different magnetic geometry. In our case all magnetic mo- 



ments and molecular field were in plane of the structure, 
whereas in their report magnetic moments were parallel 
while magnetic field was tilted out of the plane. 

We have found several interesting features in the an- 
gular and bias dependence of transport characteristics. In 
the case of an empty dot, TMR was found to be roughly 
independent of the bias voltage, but strongly dependent 
on the angle between magnetic moments. When magnetic 
moments of both leads rotate in opposite directions, both 
electric current and TMR vary nonmonotonously with in- 
creasing angle between the magnetic moments of the leads, 
and maximum of TMR may occur at a noncollinear con- 
figuration. For a singly occupied dot and for the case (i) 
(95L = 0), we found strong asymmetry in electric current 
and TMR with respect to the bias reversal, which disap- 
pears for the case (ii) (<pn = —<Ph)- This diode-like behav- 
ior of the current-voltage characteristics may be of some 
interest from the application point of view. Moreover, the 
asymmetry in current-voltage curves leads to associated 
asymmetry in TMR which may become even negative for 
one bias polarization. An important and interesting result 
is also an enhancement of TMR due to the dot level split- 
ting. Finally, we have also demonstrated numerically and 
analytically the evolution of the cotunncling gap with the 
splitting of the dot level and magnetic configuration of the 
system. 

The work was supported by the Polish State Commit- 
tee for Scientific Research through the projects PBZ/KBN/044/P03/20C 
and 2 P03B 116 25. The authors acknowledge discussions 
with Jiirgen Konig and Jan Martinek. 



A Calculation of cotunneling rate 

Here we present some details of the calculation of cotun- 
neling rates. As an example we consider the rate given by 
Eq. (4). Using the properties of the delta-Dirac function 
it can be written as 



'LR.,0 — 



L R. 



def(e)[l - /(e + ml - Mr)] 



cos 2 (<^l/2)cos 2 (vr/2) 



(e + Ml 
sin <£>l sin ip-& 

sin 2 (<^ L /2)sin 2 (^r/2)" 



Ml - £4. 



(25) 



(s + Ml -£i) 2 
Then, taking into account the identity 

f(e)[l - f(e + a = h{-Z)W + O - f(e)] , (26) 

Eq. (Al) can be rewritten as a sum of integrals of the 
type J def(s)/(e — £)" (n = 1,2,...). In principal, it is 
necessary to determine the integrals for n = 1, whereas 
the ones for n > 1 can be found using the expression 



ds 



fie) 



71+1 



1 



ds 



(27) 
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Below we demonstrate the calculation of one of the sum- 
mands, namely 



J 



de- 



e + /iL - £t 



(28) 



To calculate this integral we use the Lorentzian cutoff of 
the form g a (e) = W 2 /[(e - fi a ) 2 + W 2 }, with W being the 
cutoff parameter. Thus, Eq. (A3) can be expressed as 



J = — [ de- 

2i J (e + hl- e T )(e 



Ml - iW) 



_W r f{e) 

2i J (e + fi L - £ T )(e - fi L + iW) 



(29) 



The first (second) integral in the above formula has poles 
at £ = £f — /jl, £ = /iL + (— )iW, and e — i(2m + l)ir, with 

m = 0, 1, 2, However, because we are interested in the 

deep Coulomb blockade regime (where the second-order 
processes dominate), it is justifiable to assume e <C e a 
and this way neglect the contribution of the first pole. 
Then, by means of the contour integration and assuming 
W to be the largest energy scale, one gets 



Ml 



2 2nk B T 



hi 



W 



2nk B T 



(30) 



As the single integral depends on the cutoff parameter, 
the total rate does not. The expressions depending on W 
cancel in pairs, which can be simply seen from Eq. (A2). 
The other summands of Eq. (Al) can be found in a similar 
way with the aid of the above mentioned identities. 

Another way to calculate the cotunneling rates is to 
make use of the assumption e <C e a and approximate the 
resolvents of Eq. (Al) by l/(e + /il — £|) ~ — £t)- 

As a consequence, one arrives at the expressions of type 
J def(e)[l — f(e + £')], which can be easily calculated 
|25| . The latter method is equivalent to the former one 
if one expands the digamma functions and neglects the 
higher-order corrections in x/y, with x = |eVj,A;BT and 
y = e a ,e a + U. The advantage of using the first way of 
calculating the rates is that the higher-order corrections 
in temperature are properly described. 
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